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ABSTRACT 

B-band light-curve rise times for eight unusually well-observed nearby Type la supernovae (SNe) are 
fitted by a newly developed template-building algorithm, using light-curve functions that are smooth, 
flexible, and free of potential bias from externally derived templates and other prior assumptions. 
From the available literature, photometric BVRI data collected over many months, including the 
earliest points, are reconciled, combined, and fitted to a unique time of explosion for each SN. On 
average, after they are corrected for light-curve decline rate, three SNe rise in 18.81 ±0.36 days, while 
five SNe rise in 16.64 ± 0.21 days. If all eight SNe are sampled from a single parent population (a 
hypothesis not favored by statistical tests), the rms intrinsic scatter of the decline-rate-corrected SN 
rise time is 0.96t°;25 days - a first measurement of this dispersion. The corresponding global mean 
rise time is 17.44 ± 0.39 days, where the uncertainty is dominated by intrinsic variance. This value is 
w2 days shorter than two published averages that nominally are twice as precise, though also based 
on small samples. When comparing high-z to low-z SN luminosities for determining cosmological 
parameters, bias can be introduced by use of a light-curve template with an unrealistic rise time. If 
the period over which light curves are sampled depends on z in a manner typical of current search 
and measurement strategies, a two-day discrepancy in template rise time can bias the luminosity 
comparison by ~0.03 magnitudes. 

Subject headings: supernovas: general — cosmology: observations — distance scale 
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1. INTRODUCTION 

Within a few minutes of explosion, Type la supernovas 
release most of their energy, but due to self-absorption 
they reach peak luminosity only after 2-3 weeks. During 
this period of ballistic expansion, while the photosphere 
grows in radius but shrinks in characteristic velocity as 
slower, heavier ejecta are revealed, basic properties of 
the explosion become evident. Spectroscopic signatures 
of elements intermediate between carbon-oxygen fuel and 
iron-group ash (Filippenko 1997; Branch et al. 2006) re- 
veal that burning is incomplete, with deflagration likely 
playing an early role (Mazzali et al. 2007); nonvanishing 
polarization measures the progenitor's asphericity (Wang 
et al. 2007). After peak brightness, when iron features 
blanket the spectrum and polarizations wane, SNe be- 
come more homogeneous. For SN science to progress, 
therefore, it is crucial to study the period of rising lumi- 
nosity. 

In cosmological studies, SNe la are prized for their use 
as standardizable candles to trace the history of cosmic 
expansion (Riess et al. 1998, 2007; Perlmutter et al. 1999; 
for a review, see Perlmutter & Schmidt 2003); in this 
context, periods of greater SN uniformity are of greater 
value. Indeed, post-maximum luminosity indicators do 
yield low-dispersion Hubble diagrams (Wang et al. 2003; 
Wang, X. et al. 2005); post-maximum color measure- 
ments do solidify the corrections made for absorption by 
host-galactic dust (Lira 1995; Phillips et al. 1999; Jha 
et al. 2007). Nevertheless, as ever more ambitious cam- 
paigns to chronicle the Universe's expansion history are 
planned (e.g. Aldering et al. 2002), the fundamental issue 
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of high z — > low z evolution (Howell et al. 2007) keeps 
SN science in focus. For example, more than one SN 
la progenitor or explosion mechanism might be at work, 
with progeny neither equally bright nor equally abundant 
at high vs. low redshift. To secure such understanding, 
continued study of the rise-time period is essential. 

The subject of this report is a basic property of this 
period - the light-curve rise time itself. The i?-band rise 
time is quite sensitive to the main-sequence mass of the 
white dwarf progenitor and to its carbon/oxygen ratio 
(see e.g. Domfnguez et al. 2001). It is less sensitive to 
the progenitor's metallicity. 

Pskovskii (1984) published the first rise times for 
classes of type I SNe. In retrospect their range is reason- 
able, but, oddly, the reported correlation of rise time with 
decline rate was positive. For individual SNe, the earliest 
rise-time measurements were made by Leibundgut et al. 
(1991) (SN 1990N) and by Vacca & Leibundgut (1996) 
(SN 1994D). For a group of SNe, the earliest measure- 
ment of the average decline-rate-corrected rise time was 
reported by Groom (1998) and Goldhaber (1998). They 
used a model described e.g. by Arnett (1982), in which 
the initial rise of B flux with time is parabolic. Within ± 
w2 days, these early determinations agree with current 
values. Soon thereafter, in a definitive paper, Riess et 
al. (1999b) (henceforth Rie99b) established the presently 
accepted average rise time to B maximum of 19.5 ± 0.2 
days. This was accomplished by ferreting out early unfil- 
tered photometry for ten nearby SNe and transforming 
it to standard passbands, to which the parabolic model 
was applied. Also using this model, Conley et al. (2006) 
(henceforth Con06) recently confirmed the Rie99b low- 
redshift average, in addition measuring an average rise 
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time at high redshift of lQ.lO^J^stat) ±0.2(syst) days. 
Summarized later in this paper are additional measure- 
ments that quote larger uncertainties and/or apply to 
only 1-3 SNe; these best-fit rise times tend to be 1-3 
days shorter. 

In this paper, using literature SNe, a first analysis is 
made of the intrinsic variation of SN rise times, after cor- 
recting for their correlation with light-curve width. En- 
countered are two curiosities: a remarkable distribution 
of rise-time variations and a surprising departure from 
the presently accepted average SN rise time. Applied 
to this study is a newly developed light-curve-template- 
building algorithm that matches smooth, flexible func- 
tions to photometric data over many months including 
the earliest useful points. Externally derived templates 
and other prior assumptions are avoided. The method 
combines and reconciles the available B VRI photometry 
and fits it to a unique time of explosion for each SN. In- 
cluded in the analysis are three SNe with newly available 
early photometry. 

The bulk of this report is organized into five sections. 
In §2 the rise time is defined, including its extrapola- 
tion to the time of explosion. Described thereafter is the 
light-curve-template builder and its application to rise- 
time measurement. In §3 the SN sample and its correc- 
tions are introduced. Next a measure of SN light-curve 
width is chosen and its fitted values are discussed. The 
fitted rise times then are presented and their uncertain- 
ties are analyzed. In §4 the rise times are corrected for 
light-curve width, and their intrinsic variation is derived 
and characterized. Particular attention is devoted to the 
possibility that the parent rise-time distribution is bi- 
modal. The discussion in §5 begins by addressing the 
issue of consistency with previous average rise-time mea- 
surements. Thereafter the light-curve-width-corrected 
rise time is identified as a second SN la parameter, and 
its correlations with other variables of interest are sur- 
veyed. Finally, §6 explores the implications for precision 
cosmology of both the intrinsic scatter in the rise time 
and the present inconsistencies in its measured average. 

2. RISE-TIME DEFINITION AND FITTING 
METHOD 

2.1. Approach to Rise- Time Fitting 

Following Vacca & Leibundgut (1996), data including 
the earliest available observations are fitted to a single 
light curve in each band that everywhere satisfies strin- 
gent requirements on continuity, smoothness and analyt- 
icity. Not pursued here is the piecewise method intro- 
duced by Rie99b, in which data collected earlier or later 
than a particular cutoff time ijoin, typically 10 days be- 
fore the time £_B max of maximum B flux, are fitted to 
dissimilar functions in the two regions. In one variant 
of that method, at ij j n the two fitted curves were not 
required to meet (Rie99b; Con06). In another variant 
(Aldering et al. 2000), the two curves did meet at tj - ln , 
but the derivatives there were different. 

Most of the analysis discussed by Rie99b, and all of 
the analysis by Aldering et al. (2000) and Con06, is de- 
voted to combined data samples from ensembles of SNe, 
many of which were not measured completely enough to 
provide well-defined individual BVRI light curves with- 
out imposing external templates. Conversely, in order to 



measure the distribution of SN rise times, this analysis 
uses no templates and fits only individual SNe that do 
meet this standard. 

2.2. Characterizing Light- Curve Width 

Goldhaber (1998), Rie99b, and Con06 have established 
that SN la rise times are correlated positively with light- 
curve widths, fitting the latter using data collected no 
earlier than «10 days before t_B max (see also Riess et 
al. 1999a). In broad use are three "first parameters" 
that characterize light-curve width: Amis (Phillips 1993; 
Phillips et al. 1999); the MLCS parameter A (Riess et al. 
1995, 1996; Jha et al. 2007); and time-axis stretch sb 
(Perlmutter et al. 1997; Goldhaber, Groom et al. 2001); 
see also the parameters i±i/2 defined by Contardo et 
al. (2000). For cosmology, light-curve width is used pri- 
marily to apply the empirical broader-brighter correction 
that relates the observed luminosity of a SN to its lumi- 
nosity distance. If that correction is based on both the 
rising and falling parts of the light curve, as is true for 
A and sb, often its reliability is enhanced: measuring 
the time of peak flux, as is required to determine Amis 
without the aid of a template, is a delicate task (§3.2). 

Nevertheless, to study the correlation between light- 
curve rise time and width, it is advantageous here to 
adopt a width parameter that, like Amis, is based only 
on the declining part of the curve. Otherwise, part of the 
measured correlation would need to be attributed to the 
mutual redundancy of rise time and width. Another ad- 
vantage of Amis is that it is a simple property of the light 
curve (the gain in B magnitude between ts ma x and 15 
rest- frame days thereafter), which, for well-sampled SNe, 
can be determined without using a template (though usu- 
ally it is not). Conversely, for example, sb is measured by 
fitting the B photometry to a single rigid template (Lei- 
bundgut 1988; Goldhaber, Groom et al. 2001). When the 
photometry is exceedingly precise, the natural diversity 
of SN light-curve shapes can prevent template fits from 
achieving an acceptably probable x 2 . The best-fitted A 
or sb then may vary systematically, depending on when 
and how accurately the SN flux was sampled. For this 
study, therefore, Amis is a natural choice for measure- 
ment of SN light-curve width. 

In this paper it is convenient to rcexpress the decline 
rate Ami 5 in terms of a tightly coupled quantity that 
has the same dimensions and order of magnitude as the 
rise time t r . A natural choice is the functional inverse 
of Am.15: the fall time tf is defined as the interval after 
^Bmax (in rest-frame days) that is required for the B mag- 
nitude to dim by 1.1 magnitudes (1™1), the fiducial value 
of Amis used by Phillips et al. (1999) and many others. 
Correspondingly, the fiducial fall time is 15 days. The 
fall time is the same as the quantity i+i/2 introduced by 
Contardo et al. (2000), except that it is referenced to a 
slightly later epoch (in their notation, tf would be called 
£+0.36)- Over the range of Ami 5 appropriate to SNe dis- 
cussed in this paper, tf is > 99.9% correlated with Amis; 
the two variables are related by 

t f « 15-8.42(Amis-l.l) + 4.54(Amis-l.l) 2 . (1) 

Also used here is the rise-time - fall-time difference t T { = 
t r - tf. 

Occasionally, for displaying and discussing the prop- 
erties of B and V light curves, this paper uses a time 
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dilation factor s = tf/15 days that likewise is > 99.9% 
correlated with Amis. This factor is used only to define 
the dilated rest-frame phase r, 

so that B(t=15 days) - B(t = 0) = l m l. Conveniently, 
when B is plotted vs. r, all light curves exhibit the fidu- 
cial decline rate. 

Finally, to make a rough estimate of the sensitivity of 
existing cosmological analyses to uncertainties in the rise 
time of the light-curve templates they have adopted, it is 
proposed (§6.2) that the stretch sb that is measured in 
these analyses is approximately proportional to a linear 
combination t r + jtf of rise and fall times. That is, for an 
ensemble of SNe, a plot of t T or if vs. sb would scatter 
about a straight line with finite slope. This stands to 
reason: all three variables are amplified by the same fac- 
tor if the time axis is dilated. Because the declining part 
of the light curve is longer and may have smaller photo- 
metric errors than the rising part, one might expect sb 
to be more strongly correlated with tf than with t T ; cor- 
respondingly, 7 is assigned the value 2 in the examples 
given. 

2.3. Extrapolating to Time of Explosion 

Conventionally, the SN la rise time f r is defined as the 
interval between the time of explosion t exp \, when fluxes 
are assumed to vanish, and i_B ma x- When this defini- 
tion is adopted, as it is here, measuring t r requires ex- 
trapolating to t eX pi- This is necessary because available 
template-quality filtered photometry for nearby SNe be- 
gins at most «3 magnitudes below peak. If the initial 
time dependence of the B flux is parabolic (Arnett 1982), 
for typical B light curves t CKp i occurs at least ~3 days 
before the first datum. 

When a particular function such as a parabola is cho- 
sen for this >3-day extrapolation in the B band, the 
corresponding function for other bands can have exactly 
the same form only if the SN color does not evolve within 
the region of extrapolation. Unfortunately, as is exhib- 
ited clearly in Fig. 7(b) of Pastorello et al. (2007b), this 
assumption is not fully supported by data available at 
slightly later epochs. To explore this point, consider the 
rate of change vx = d(B — X)/dr of B — X color, where 
X = V, R, or /, in the period r < — 11 days. In this 
analysis, vx is parametrized by a constant, whose value 
is determined primarily by smooth extrapolation of the 
fit to B — X data in the region —11 < r < — 5 days. 
Considering as diverse examples SN 1998aq, SN 1994D, 
SN 2002bo, and SN 2005cf, respectively, the best-fit v v 
is approximately m 000, -0"?011, -0 m 045, and -0 m 067 
dy _1 . For the last three SNe, the departure from vv = 
is significant and nonnegligible compared to the maxi- 
mum value d(B— V)/dr w +0 n ?06 dy" 1 typically reached 
during the post-i3 max period 15 < r < 20 days. If the 
B-band flux rises parabolically, 

f B oc (r - r cxpl ) 2 , (3) 

the flux in the X band will exhibit a slightly slower rise, 

/xoc(r-r cxpl ) 2 10 - 4 ^^-^, (4) 

when vx is negative. 



Using a parabola for the B-band extrapolation to t e xpi 
is only one possible method for deriving t r from the ob- 
served data. It is a simple and conventional method, ap- 
plied uniformly to each SN, that conforms to the earliest 
available data. Should it be discovered that, at epochs 
too early to be well measured at present, the B flux does 
deviate from parabolic time dependence, rise times that 
are corrected for that effect would shift systematically 
relative to those reported here. 

2.4. The AQUAA Template Builder 

For constructing light-curve templates, a primary chal- 
lenge is to devise functions that approach the smooth- 
ness of the parametrization used by Vacca & Leibundgut 
(1996), while attaining the flexibility required to follow 
the precise photometry in detail. To achieve this flexibil- 
ity, in earlier approaches cubic splines have been fitted to 
magnitudes (Hamuy et al. 1996; Prieto et al. 2006) or to 
fluxes (Knop et al. 2003). Because magnitudes soar near 
the time of explosion, the former type of fit has found 
use mainly after is ma x ~ 5 days; because splines do not 
naturally mimic exponential decay, the latter type of fit 
has proved awkward in the region beginning «50 days 
after i_B max , requiring an exponential to be spliced in. 
As well, if the cubic spline knots are placed arbitrarily, 
unwanted ringing can arise. 

Since August 2005 the template builder AQUAA 
(Adaptive QUArtic Algorithm) has been under develop- 
ment, with the primary aim of producing full (U)BVRI 
light-curve templates representing individual SNe (Stro- 
vink, M. 2007, in preparation). AQUAA is based on fits 
to quartic splines, which are smoother than cubic splines 
(third derivatives are continuous). Further smoothing is 
gained by using "fuzzy knots" that allow gradual tran- 
sitions between spline segments; the resulting curves 
are all-orders differcntiablc. Knot positions are deter- 
mined adaptively and objectively by the fit itself, sub- 
ject only to minimum separation requirements that are 
SN- independent. 

A quartic spline of this type is used to model the mea- 
sured i?-band flux. In its earliest segment, a parabolic 
rise is imposed by requiring other coefficients to vanish. 
Near-exponential decay in the region «50 < r < 125 
days is modeled by using as the independent variable 

u = m tanh ((t - t B max)/w ) , (5) 

where uo « 125 days, in place of the time t itself. In the 
spline's latest segment, the cubic and quartic coefficients 
are set to zero. 

To fit the VRI bands, 1 similar quartic splines are used 
to model the B — V, B—R, and B—I colors. This choice is 
made because color curves have finite asymptotes; they 
vary less violently than magnitude curves; and, in the 
/ band, they have fewer extrema. Implicitly it is as- 
sumed that colors evolve smoothly even near the (com- 
mon) time of explosion, where magnitudes change vio- 
lently; when early filtered data are scarce, this assump- 
tion enhances the fitted rise-time accuracy. All coeffi- 
cients for all splines are fitted simultaneously by mini- 
mizing a global x 2 that compares modeled to measured 

1 For a few SNe with sufficient data, the U—B color may similarly 
be modeled. The U band is not used here for rise-time determina- 
tion. 
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B VRI magnitudes and their quoted errors (assumed in- 
dependent band-to-band) . Because the spline coefficients 
are overconstrained by the asymptotic requirements, the 
global fit requires the use of Lagrange undetermined mul- 
tipliers. 

Typically an average of 9 degrees of freedom (dof) per 
band are determined by the global fit. In advance, scores 
of simpler fits are performed in order to zero in on the 
global optimum. Currently AQUAA requires at least one 
day of human effort and CPU time to fit one SN. 

2.5. Applying AQUAA to Rise-Time Measurement 

Figure 1 shows AQUAA fits to the eight literature SNe 
fitted for this rise-time study. Although the rise time 
is the product of a global fit to all four BVRI bands 
over available data out to r = 125 days, Fig. 1 exhibits 
only B-band fits to SN 2005cf, SN 2003du, SN 2002bo, 
SN 1994D, SN 2001el, and SN 2004eo - and only V fits 
to SN 1998aq and SN 1990N - over r < 17 days. For 
the latter two SNe, the V band is shown in order to 
include the unfiltered points transformed to that band 
by Rie99b; the fitted U-band maxima are reached 1.2 
and 1.9 rest-frame days, respectively, after iB max - On 
the ordinate, the square root of the flux (normalized to 
its t = value) is displayed so that a parabolic rise 
from explosion appears as a straight line. For ease of 
comparison, all SNe are forced to have Ara r= i5 = 1™1 
by using r from equation (2) as the abscissa. Therefore, 
by construction, all B light curves agree at r — and t = 
15 days. Nevertheless, close inspection reveals detailed 
differences in the shapes of the declining parts of the 
curves: without constraint from any template, AQUAA 
draws an independent smooth curve through each SN's 
data. 

Another feature evident from Fig. 1 is that all fits 
are "good" (x 2 ~ #dof). This also occurs in part by 
construction. For each combination i of telescope and 
band that contributes to the SN dataset (telescopes are 
grouped if many contribute), both the zero point and 
the photometry error may be perturbed. To perturb the 
zero point, an offset Si and an offset tolerance are in- 
troduced; initially all Si = and all are equal. The 
AQUAA fitted x 2 i s augmented by ^2A$i/ei) 2 and the 
optimum {Si} are determined as part of the global fit. 
This process is repeated several times, with the {ei} ad- 
justed iteratively so that the offset contribution to \ 2 
is balanced among telescopes and bands, with a total 
of « l/#dof. When only two telescopes are involved, 
this procedure assigns their zero-point | offsets | equally; 
but when many are involved, the iteration assigns tighter 
tolerances to telescopes that are in better mutual agree- 
ment. Typically, telescopes with significant influence on 
the fitted light-curve shapes are assigned offset tolerances 
of 0™02 or less. Rarely, when the fit strongly prefers it, 
a small color-term offset is combined with the zero-point 
offset. To perturb the photometry error, an error floor 
(typically of order 0™015) may be imposed, and/or a con- 
stant error (typically of the same order) may be added in 
quadrature to the quoted uncertainties. Again the pro- 
cedure is iterated, assigning the largest error augmenta- 
tions to the telescopes with greatest scatter (relative to 
quoted errors), while leaving unmodified the errors on 
the smoothest data. 

As for the rising part of the light curves, on the scale 



of Fig. 1 the root fluxes in the B band grow in a manner 
that is indistinguishable from linear over the first «40% 
of their rise times. Thereafter, everywhere maintaining 
a continuous third derivative, the curves gradually turn 
over in order to match the data around the peak. Close 
inspection reveals, however, that the U-band root flux 
fitted to SN 1990N has a slower-than-linear rise. As was 
discussed in §2.3, this occurs because the initial rate of 
change vy oi B — V color, based on the best AQUAA fit to 
the nearest available SN 1990N data, is significantly less 
than zero. As will be discussed in §3.4, uncertainties in 
the the AQUAA fitted output take into account the error 
associated with this extrapolation. Appendix A supplies 
more detail on fits carried out for individual SNe. 

The striking feature of Fig. 1 is that the fitted values 
of t at which the explosions occurred, indicated by "x", 
appear to divide into two groups separated by ~2 days. 
The diversity of decline-rate-corrected rise times would 
be apparent even if no fitted curves were displayed; to 
this topic the discussion now turns. 

3. RISE TIMES OF EIGHT TEMPLATE SNE 
3.1. SN Sample Characteristics 

At present, individual SNe with published photometry 
having sufficient cadence and accuracy to define indepen- 
dent BVRI templates are available only at low z (< 0.02). 
The eight nearby literature SNe studied here are listed 
in the first column of Table 1. Five of these (SN 1990N, 
SN 1994D, SN 1998aq, SN 2001el, and SN 2002bo) are 
among the eight nearby SNe whose rise times were ana- 
lyzed by Con06. The first three in this list are among the 
ten SNe studied by Rie99b; in particular, SN 1990N had 
a large influence on that paper's rise-time results. The 
analysis reported here docs not use the four SNe studied 
by Rie99b but not by Con06, nor docs it use three of 
the SNe studied in both papers. For six of these seven 
SNe, the quality of the published photometry does not 
allow AQUAA to define BVRI templates independently 
for each SN. For the seventh, SN 1998bu, the only pub- 
lished datum earlier than «4 days before t s ma x is a single 
unfiltered point with a sizable uncertainty (~ 0™15). 

To the five SNe in common with Con06, this analysis 
adds three that are more recently measured: SN 2003du 
(Leonard et al. 2005; Anupama et al. 2005; Stanishev 
et al. 2007; Li, W.D. et al. 2007, in preparation); SN 
2004eo (Hamuy et al. 2006; Pastorello et al. 2007a); and 
SN 2005cf (Pastorello et al. 2007b). These recent SNe 
are splendidly observed and contribute greatly to this 
report. 

Excluding SN 2004eo, this sample has an average he- 
liocentric redshift (z) = 0.004 with an rms a z — 0.002; 
such a low g z introduces a negligible rise-time dispersion. 
Following usual practice (see, for example, Prieto et al. 
(2006)), SNe with z < 0.01 were not if-corrccted; SN 
2004eo, with z = 0.016, was ^-corrected. To first or- 
der, to compensate for dust absorption Ax in band X, 
dereddening (R-) corrections add a constant offset — Ax, 
which has no effect on this analysis. To second order, 
i?-corrcctions perturb the light-curve shape. This oc- 
curs mainly where colors are changing rapidly, affecting 
the B light curve's falling part more than its rising part. 
SN 2002bo and SN 2001el were i?-corrcctcd using the 
A(X) of Cardelli et al. (1989). The other six SNe, having 
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Fig. 1. — Representative AQUAA light-curve fits. The abscissa r is the SN phase relative to tB max , dilated by a factor s so that 
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((B-V)(t Bmax )) = -0 m 01 with an rms of m 08, on av- 
erage are not too far from (B — V)(tBmax) ~ — 0™05, as 
is typical of fiducial SNe la, and so were not i?-corrected. 
The B and V data from Fig. 5 of Hamuy et al. (2006) 
that contribute to the AQUAA fit of SN 2004eo were (S-) 
corrected to the standard filters of Bessell (1990) using 
the procedure of Stritzinger et al. (2002). 

In summary, three of the eight SNe received a correc- 
tion: SN 2001el (R); SN 2002bo (R); and SN 2004eo 
(KS). The combined (KRS-) corrections were calcu- 
lated using vl.2 of the spectral template series introduced 
by Nugent et al. (2002). 2 Iteratively these templates were 
dilated and warped so that, after X.RS'-corrections, the 
synthesized B — V color evolution matched the observed 
evolution. Within the range of Fig. 1, the largest \KRS- 
corrections| to B-B max for SN 2001el, SN 2002bo, and 
SN 2004eo, respectively, were small: -0"?044, -0™064, 
and (owing to a cancellation) +0™009. 

All SNe that ever have been fitted by AQUAA are in- 
cluded in the rise-time sample, except for SN 1995D, 
SN 1996X, and SN 2002fk. The first two are not use- 
ful for rise-time analysis due to lack of _B-band data for 
t < — 4 days. For the last, which is similar in rise and 
fall time to two SNe in the rise-time sample, the pho- 
tometry is not yet published (Li, W.D. et al. 2007, in 
preparation). Therefore the sample is unbiased by fore- 
knowledge of fitted rise-time values. SNe belonging to 
the subsets "1991T-like" and "1991bg-like" (Filippenko 
et al. 1992a, b) that are equally well suited for rise time 
analysis were not identified and therefore not added to 
the sample. Excluded were more peculiar SNe such as 
SN 2000cx, which exhibited unusual color evolution and 
spectral features (Li et al. 2001) along with an obviously 
short rise time. The result is a standard SN la set similar 
to one that might be used in a precise cosmological study. 
Its members are unusually well observed, possessing nor- 
mal characteristics and requiring very little correction. 

3.2. Fitted Decline Rates and Fall Times 

Care must be exercised in fitting the decline rate Amis: 
because the time derivative of the flux at day 15 is large, 
a small bias in the definition of ts ma x can propagate into 
a large bias in Ami 5 . The AQUAA fitted decline rates 
appear in column 3 of Table 1. Whenever Amis values 
and uncertainties are quoted in the photometry refer- 
ences tabulated there, they are found to be in accept- 
able agreement with the AQUAA fitted values. For both 
SN 1990N and SN 1994D, similar concordance is found 
with values quoted both by Phillips et al. (1999) and by 
Rie99b. 

In Fig. 2 the filled squares show the AQUAA fitted fall 
times tf vs. ATO15 for this sample (with SN 1995D and 
SN 1996X added). The exhibited best-fit quadratic curve 
through these points is given by equation 1 (§2.2). When 
the ordinate is transformed so that this curve becomes a 
straight line, the filled squares in the transformed Fig. 2 
fully correlate the ordinate with the abscissa (Pearson 
R > 0.999). 

3.3. Fitted Rise Times 

The AQUAA fitted rise times appear in column 4 of Ta- 
ble 1 and as filled circles in Fig. 2. In addition, plotted 

2 http : //supernova. lbl . gov/~nugent 



as an open triangle is the average fitted rise time for ten 
SNe corrected to unit stretch by Rie99b; shown as an 
open square is a similar average for the eight low-z SNe 
studied by Con06. These two points are plotted at the 
fiducial Amis values that corresponded to unit stretch 
for each analysis (Rie99b; Conley, A., private commu- 
nication). Both statistical and systematic uncertainties 
are included in all error bars (see §3.4). 

As is discussed in Appendix B, for five SNe the rise 
times in Table 1 may be compared to published deter- 
minations that used the piecewise method of Rie99b (see 
§2.1). Compared to AQUAA values, the piecewise method 
yields rise times that, on average, are longer by 1.28±0.33 
days. Of this systematic difference, 0.44 ± 0.15 days are 
due to the later ts ma x that was used by the published 
rise-time determinations. The extreme instance of this 
(ifimax later by 0.9 days for SN 1990N) may be ascribed 
to the use by Rie99b of only a second-order polynomial to 
fit the B peak. If a higher-order polynomial is used, the 
empirical fact that B light curves rise faster than they 
fall causes the fitted peak to become slightly asymmet- 
ric, yielding a slightly earlier i_B ma x- The balance of the 
systematically longer rise time obtained by the piecewise 
method occurs because that method assumes that B — V 
color does not vary within the early piece, and because 
it does not require the two pieces to join. This point is 
elaborated in Appendix B. 

As foreshadowed by Fig. 1, the striking conclusion from 
the AQUAA fitted rise times and decline rates plotted in 
Fig. 2 is that SN la rise times diverge significantly, in 
a way for which one cannot compensate by applying a 
purely decline-rate-dependent correction. 

3.4. Error Analysis 

A distinct advantage of using a single global fit to all 
the relevant (BVRT) data to determine a parameter such 
as the rise time is the simplicity of the (frcquentist) ap- 
proach used here to propagate photometric errors into 
the fitted parameter's uncertainty. In the global fit, when 
the parameter in question is varied and the change in like- 
lihood is recorded, all other (m 4 x 9) free parameters are 
allowed also to vary for best compatibility with that ex- 
cursion. (One such free parameter, defined in §2.3, is uy, 
the coefficient of the linear term in the earliest segment 
of the quartic spline representing B — V color.) Then 
the curvature at maximum of the one-dimensional log 
likelihood function yields the photometric error on the 
parameter in question. In this approach, as for all anal- 
ysis in this report, no Bayesian priors are employed. (If, 
alternatively, an integral over a multidimensional likeli- 
hood surface had been performed, implicitly a Bayesian 
prior would have needed to be applied.) 

As is often the case, here the most challenging aspect 
of error analysis is not the propagation of random errors, 
but rather the estimation of systematic errors. In the 
development of AQUAA, many algorithmic choices were 
made. For example, a variety of schemes for requiring 
a minimum degree of light-curve smoothness were ex- 
plored. Several of these schemes are plausible and de- 
fensible; nominally, choosing among them should have 
a neutral effect upon any fitted parameter. Eventu- 
ally a particular scheme was adopted. However, despite 
their nominal neutrality, such algorithmic choices had 
discernible impacts on some fitted values. Fortunately, 
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TABLE 1 

Results of AQUAA light-curve fits. 



SN 




^Bmax 
(JD - 2,440,000) 


Amis 
(mag) 


Rise time 
(rest-frame 
days) 


Rise time 
— fall time 
(rest-frame 
days) 


Photometry 
references 


SN 1990N 




8082.46 ± 0.35 


0.990 ±0.034 


20.01 ±0.46 


4.04 ±0.74 


(1,2) 


SN 1994D 




9432.29 ± 0.14 


1.344 ±0.021 


15.39 ±0.47 


2.04 ±0.53 


(3,4) 


SN 1998aq 




10931.04 ± 0.14 


1.042 ±0.021 


17.52 ±0.58 


2.03 ± 0.64 


(5,2) 


SN 2001el 




12182.19 ± 0.17 


1.168 ±0.021 


18.00 ±0.56 


3.57 ±0.63 


(6,7) 


SN 2002bo 




12356.39 ± 0.20 


1.162 ±0.031 


16.05 ±0.37 


1.57 ±0.53 


(8,9, 10) 


SN 2003du 




12765.95 ± 0.14 


0.982 ±0.021 


17.71 ±0.35 


1.64 ±0.44 


(11, 12, 13, 14) 


SN 2004eo 




13278.65 ± 0.14 


1.403 ±0.037 


16.64 ± 0.44 


3.86 ±0.53 


(15, 16) 


SN 2005cf 




13533.78 ± 0.14 


1.068 ±0.021 


16.62 ±0.25 


1.35 ±0.36 


(17) 


Fiducial (slow rising a ) 




1.1 


18.81 ±0.36 


3.81 ±0.36 


Fiducial (fast rising b ) 




1.1 


16.64 ±0.21 


1.64 ±0.21 




Fiducial (all c ) 






1.1 


17.44 ±0.39 


2.44 ± 0.39 





References. — (1) Lira et al. (1998); (2) Riess et al. (1999b); (3) Richmond et al. (1995); (4) Patat et al. 
(1996); (5) Riess et al. (2005); (6) Krisciunas et al. (2003); (7) Krisciunas et al. (2007); (8) Benetti et al. (2004); 
(9) Krisciunas et al. (2004); (10) Szabo et al. (2003); (11) Leonard et al. (2005); (12) Anupama et al. (2005); (13) 
Stanishev et al. (2007); (14) Li, W.D. et al. (2007), in preparation; (15) Hamuy et al. (2006); (16) Pastorello et 
al. (2007a); (17) Pastorello et al. (2007b). 

a Average of SN 1990N, SN 2001el, and SN 2004eo, corrected to Amis = l m l. 
b Average of all but SN 1990N, SN 2001el, and SN 2004eo, corrected to Amis = 

c Average of all SNe, corrected to Amis = 1™1, and with errors augmented to account for intrinsic variance. 



over the long AQUAA development period, these fitted 
values were recorded regularly, for example since April 
2006 for SN 1994D and SN 2001el. For most SNe, these 
records make it possible, for example, to estimate a sys- 
tematic rms ismax and rms Amis. Conservatively, for 
a particular parameter, each SN is assigned a system- 
atic error equal to the median systematic rms for that 
parameter, or its own, whichever is greater. 

Particular attention has been devoted to systematic er- 
ror on the rise time t T . First, it should be emphasized 
that any systematic uncertainty in i_B ma x propagates di- 
rectly into a systematic error in t v (and in fall time tf). 
If tf is subtracted from t T , as will be discussed in §4.1, 
this source of error doubles. As a further check, for all 
but SN 200 lei and SN2003du, the extent and quality of 
flux measurements at r < —10 days permit a meaningful 
parabolic fit to those data alone. This cross-check was 
carried out in the B filter only, except for SN 1990N and 
SN 1998aq, for which the F-filter data that largely con- 
trol their times of explosion i cxp i were included. From 
these simple fits, t cxp \ on average differed by only +0.01 
days from those of AQUAA, with an rms difference of 
0.33 days. Conservatively, for these SNe, the full dif- 
ference in texpi between the simple and AQUAA fits was 
added in quadrature to the other rise-time errors; for 
the remaining SNe, the median of those |differences| was 
used. Finally, the possibility of Malmquist bias in the 
measurement of t r was discussed by Rie99b. As will be 
seen, the analysis here is focused mainly on the B light 
curve's symmetry, which to leading order is luminosity- 
independent. 

In summary, the parameter errors listed in Table 1 
and displayed in Figs. 2 and 3 are the quadrature sum 
of (i) straightforwardly propagated photometric random 
errors; (ii) systematic errors estimated from records of 
AQUAA algorithm evolution; and, for rise times, (Hi) dif- 
ferences in £ e xpi between AQUAA and simple parabolic fits 
to early data. 
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Fig. 2. — Fitted rise time or fall time (defined in the text) vs. de- 
cline rate Amis • Filled circles show the AQUAA fitted rise times for 
individual SNe. For each SN, both statistical and systematic un- 
certainties are included in its horizontal and vertical error, which 
are positively correlated. Filled squares show the AQUAA fitted fall 
times for the same SNe (plus two others with scant pre-maxirnum 
data). The quadratic curve is drawn to guide the eye. Open sym- 
bols show the results of previous rise-time fits by Rie99b (triangle) 
and Con06 (square) to ensembles of stretch-corrected low-z SNe. 

4. STATISTICAL PROPERTIES OF THE 
RISE-TIME SAMPLE 

4.1. Correcting the Rise Time for Decline Rate 

As noted in §2.2, Goldhaber (1998), Rie99b, and Con06 
independently found that SN la rise times are correlated 
positively with light-curve widths. The anticorrelation 
of rise time with decline rate that is apparent in Fig. 2 
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supports these findings (although this support would lose 
statistical significance if SN f990N were dropped as an 
outlier). In the analysis that follows (which retains SN 
1990N), the average rise time is referenced to a fiducial 
decline rate, and the portion of the rise-time variation 
that is independent of decline-rate variation is studied. 
For these purposes, the functional dependence of rise 
time on decline rate needs to be characterized. How- 
ever, it is nontrivial to fit even a simple function to 
data whose abscissae and ordinates are not highly cor- 
related, with errors that are nonncgligible fractions of 
their ranges (see, for example, the discussion by Wang, 
Strovink ct al. (2006) of the fit to the data in their Fig. 6). 

The approach taken here is to apply only a simple ad 
hoc correction for this anticorrelation - a choice that is 
not necessarily optimal, but that does, within the statis- 
tical error, obviate the need for further correction. The 
form chosen for this correction is motivated by the im- 
pression that, in Fig. 2, the rise time and fall time ex- 
hibit a similar decline with increasing A77115. Therefore, 
to correct the rise time for decline rate, the fall time is 
simply subtracted from it. This difference t r { = t r — if is 
exhibited in the fifth column of Table 1, and it is plotted 
vs. Amis m Fig. 3(a) and (b). As seen in Fig. 3(b), the 
points there do not substantially correlate t r f with A77115 
(Pearson R = 0.25±g;^). Because of this feature, the sta- 
tistical properties of f r f can be studied without reference 
to a particular decline rate: the analysis discussed below 
assigns to t T f a Amis-independent mean and dispersion. 

4.2. Intrinsic Variation 

As presented in Table 1 and shown by the error bars 
in Fig. 3(a) and by the inner error bars in Fig. 3(b), the 
total error on t T { is substantially larger than on t r . As 
discussed in §3.4, this arises in part from the doubled 
effect of the uncertainty in tsmax upon t T {, compared to 
its effect upon t r . Notwithstanding the larger error, an 
attempt to assign a unique best-fitted t T { using these er- 
rors would yield an unacceptably high x 2 (29.4/7 dof, 
confidence level 1.2 x 10~ 4 ). Therefore, beyond the total 
uncertainty cr™ cas with which it is measured, t r f demon- 
strates a nonvanishing intrinsic variation ctJ."*. The sum 
in quadrature of ctJ."* and cr^ cas is shown by the outer er- 
ror bars in Fig. 3(b), which are dominated by ctJ."*. The 
quoted value 

ajf = 0.96±g:|§ days (6) 

is obtained by varying ctJ."* to yield confidence levels of 
0.5, 0.84, and 0.16, respectively, for fits to unique values 
of f r f. 

4.3. Sample Average 

Using this central value of trj." * , the best- fitted average 
difference of rise time and fall time is 

(i rf ) = 2.44 ± 0.39 days, (7) 

where the uncertainty is dominated by the intrinsic vari- 
ance arising from crj."*. This value and its uncertainty 
are illustrated by the horizontal lines in Fig. 3(b). Cor- 
respondingly, a hypothetical large sample of SNe la ex- 
hibiting the fiducial fall time of 15 days would have a 
best-fitted average rise time 

{tf d ) = 17.44 ± 0.39 days. (8) 



4.4. One Parent Population or Two? 

In Fig. 3(a), separate values of (t r f) are fitted to the 
three slowest-rising and the five fastest-rising SNe, with- 
out augmenting the measurement errors to account for 
any intrinsic variance. Obviously these fits are satisfac- 
tory. The separate averages are 

slow: (t r f) = 3.81 ± 0.36 days; 
fast: (trf) = 1.64 ±0.21 days. (9) 

Belaboring this point, Fig. 4 displays an unbinned his- 
togram (defined in its caption) of the eight individual 
determinations of t T { . Despite the meager statistics, it is 
difficult to view Fig. 3(a) or Fig. 4 without considering 
the possibility that these SNe might be drawn from two 
different parent populations. 

The present statistics do not support an exhaustive 
attempt to analyze all the circumstances that might ac- 
count for the distributions in Figs. 3 and 4. For Monte 
Carlo simulation, only two simple models are considered 
here. In model A, 3 SNe arise from one parent popu- 
lation, and 5 from a different population. Each parent 
population has its own unique true value of t r {. The dif- 
ferences between observed and true i r f are drawn from 
normal distributions with rms values taken from the ver- 
tical error bars in Fig. 3(a). Conversely, in model B 
all 8 SNe arise from the same parent population, again 
having a unique true value of t r f. Here the differences 
between observed and true t r { are drawn from normal 
distributions using the outer error bars in Fig. 3(b). In 
the simulated analysis of each 8-SN experiment, two free 
parameters are optimized for either model. For model 
A, two values of (t T {) are fitted; for model B, one value 
of (trf) is fitted, and cr^f* 1S adjusted so that, after it is 
added in quadrature to the inner error bars in Fig. 3(b), 
X 2 = 1 /dof for that fit. 

The present SN sample cannot discriminate definitively 
between models A and B. Not only are the statistics in- 
adequate, but these models and their statistical tests are 
contaminated by having been devised after the sample's 
statistical behavior had become evident. Such a posteri- 
ori estimates of relative likelihood are notoriously unreli- 
able. Nevertheless, because of this second issue, it is still 
valuable here to discuss statistical tests of models A vs. 
B: in the future, when statistically independent samples 
of SNe for which individual rise times can be measured 
become available, the same models and tests described 
here may be reapplied without such contamination. 

Figure 5 illustrates two statistical tests of the rela- 
tive likelihood that the actual eight-SN sample is de- 
scribed by model A vs. model B. In this figure, the 
four plots on the left (right) refer to model A (B). The 
crosses in Figs. 5(a)-(b) show the distribution of the data 
in pull (the difference between the observed i r f and its 
best-fit average, divided by the total uncertainty in t T {). 
The points (with invisibly small uncertainties) show the 
Montc-Carlo-simulated pull distribution of SNe fitted in 
accord with one of the two models, and the curves are 
gaussians with the same rms width as the distribution 
of points. Because x 2 is additive, the simulated pull dis- 
tribution for model A is exactly gaussian. For model B, 
however, the ct™ * that is chosen for each simulated eight- 
SN experiment is affected by fluctuations in the values 
of trf that were generated for that experiment; therefore, 
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Fig. 3. — (a) Difference between the rise time and fall time shown in Fig. 2 vs. decline rate Amis. The error bars represent the 
combined statistical and systematic uncertainties. Horizontal lines depict the ordinates best fitted to the slowest-rising three SNe (top) 
and to the fastest-rising five SNe (bottom) with their one-standard-deviation uncertainties. The exhibited Pearson coefficients R show that 
the abscissa and ordinate are not significantly correlated for either set of points, (b) Same as (a) except that a single ordinate is fit to all 
eight SNe using the outer error bars, which include a 0.96-day intrinsic rise-time variation added in quadrature. 
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Fig. 4. — Unbinned histogram of differences t r f between rise and 
fall times for 8 SNe. Each SN is assigned a box of unit area and of 
width equal to its photometric FWHM resolution. The unbinned 
histogram is the sum of the boxes. Superimposed are expected 
summed probability density functions based on two different par- 
ent distributions convolved with the photometric resolutions. The 
parent distributions are (model A) two <5-functions in i r f , represent- 
ing 5 fast- and 3 slow-rising SNe; and (model B) a single gaussian 
of width o-jf. 

as confirmed by Fig. 5(b), the expected pull distribution 
is only approximately gaussian. 
The usual estimator of population kurtosis, 3 

G ^ ( n -"2)(n 1 - 3 ) ('" + 1) 5- 3( "-")- < 10 > 

where m 2 is the sample variance and is the fourth 
sample moment about the mean, reaches a minimum for 
distributions consisting of two ^-functions. (A sample of 
n = 8 values divided 4:4 between two <5-functions has 
G 2 = -2.8, while one divided 5:3 has G 2 = -2.24.) 
Therefore G 2 is a natural estimator of a distribution's 
bimodality. When analyzed in accord with model A, 

3 http: //en.wikipedia. org/wiki/Kurtosis 



the eight actual SNe have G 2 = —0.31 - a typical value 
- while according to model B they have G 2 = —1.98. 
These values appear as arrows in Figs. 5(c)-(d), to be 
compared to the simulated kurtosis distributions shown 
as histograms. For model B, 26 in 1000 simulated ex- 
periments yield a kurtosis at least as small as the one 
obtained for the actual sample. 

The second statistical test requires mapping the pull 
into a "pull probability" variable that is distributed uni- 
formly in the interval (0,1). These uniform distributions 
are shown for simulated SNe (points) in Figs. 5(e)-(f); 
again the eight actual SNe appear as crosses. For model 
B, between the fifth- and sixth-ranked actual SNe there 
is a large gap (0.47) in pull probability. Such a gap is 
an estimator of bimodality that is different from (but 
not independent of) the kurtosis. Figures 5(g)-(h) dis- 
play histograms of the simulated distribution of maxi- 
mum gap between adjacent pull probabilities, again with 
the actual values shown as arrows. To construct these 
maximum gap distributions, gaps between SNe ranked 2 
& 3, 3 & 4, 4 & 5, 5 & 6, and 6 & 7 in pull probability 
were considered. For model A, 957 in 1000 experiments 
would yield a larger gap; this may suggest that the sys- 
tematic errors on t T { were estimated too conservatively. 
For model B, 14 in 1000 simulated experiments yield a 
gap in pull probability at least as large as the one ob- 
tained for the actual sample. 

5. DISCUSSION 

Given the meager statistics, the topics chosen for the 
remaining discussion are relevant to any number of par- 
ent populations. 

5.1. Consistency of Measured Average Rise Times 

For SNe la with the fiducial decline rate, the AQUAA 
average rise time (tf d ) = 17.44 ± 0.39 days is m2 days 
shorter than the stretch-corrected value 19.42 ±0.19 days 
reported by Rie99b. (Two other decline-rate-correction 
methods employed by those authors yielded essentially 
the same result.) It is «2.1 days shorter than the stretch- 
corrected low- .z value 19. 58^^19 days reported by Con06, 
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Fig. 5. — Comparison to gaussian expectation of the distribution of pulls (residuals normalized to their uncertainties) for the fits shown 
in Figs. 3(a) (left column) and 3(b) (right column), (a)-(b) Distribution of pulls for eight actual SNe (crosses) and for fits to eight simulated 
SNe (points). The gaussian curve has the same rms as the points, (c)-(d) Kurtosis of the pull distribution for eight actual SNe (arrow) 
and for sets of eight simulated SNe (histogram), (e)-(f) Same as (a)-(b) except that the pull is mapped to a function ("pull probability") 
that is distributed uniformly on (0,1). (g)-(h) Distribution of the maximum interval in pull probability between neighboring interior SNe, 
for eight actual SNe (arrow) and for sets of eight simulated SNe (histogram). 



and it is ~1.7 days shorter than the high-z rise time of 
19-lOig 1? (stat) ± 0.2 (sys) days reported there. Using a 
well-observed nearby sample of similar size, and apply- 
ing these observations efficiently to the measurement of 
(t^ d ) , AQUAA nevertheless yields twice the low-z rise-time 
uncertainty quoted in those papers. 

About 2 /3 of the inconsistency in central value is due 
to the tendency, discussed in §3.3 and in Appendix B, 
for the piecewise method of Rie99b to yield rise times 
longer than those obtained by AQUAA. The balance of 
the central- value difference, as well as most of the differ- 
ence in assigned uncertainty, is due to intrinsic variance. 
Implicitly, it was assumed by Rie99b, and for low z by 
Con06, that particular small samples, in which individ- 
ual SNe such as SN 1990N carry heavy weight, accurately 
mirror the general SN la population; no error contri- 
butions from intrinsic variance were assigned. As well, 
when a particular externally derived light-curve template 
was applied by those authors, implicitly it was assumed 
to represent accurately the average SN la behavior. In 
fact, B templates in wide use differ substantially in their 
asymmetry about i Smax : B(t Bmax + At) - B(t Bm ^- At) 
can vary by up to one-half magnitude at At ps 10 days. 



Conversely, the analysis reported here applies no exter- 
nally derived template, and its uncertainty on the aver- 
age rise time is dominated by intrinsic variance. 

Among additional determinations of SN la rise time, 
at low z Aldering et al. (2000) obtained the value 
20.08 ±0.19 (stat) days. Very recently, Garg et al. (2007) 
reported a F-band rise time of 17.6±1.3 (stat)±l.l (sys) 
days based on broadband observation of 3 SNe with red- 
shifts 0.135 < z < 0.165. Since V max usually occurs 
~1.5 days later than t Bmax , their central value would 
correspond to a B-band rise time of only «16 days. At 
high z, Groom (1998) reported an average rise time of 
17.6 ± 0.4 days based on a light curve to which 37 SNe 
contributed. Applying a more complete error analysis 
to a large subset of those SNe, Aldering et al. (2000) 
obtained t r = 18.3 ± 1.2 (stat) if'J ( S Y S ) days, where the 
systematic errors are upper limits; see also (Goldhaber, 
Groom et al. 2001). When systematic as well as sta- 
tistical errors are taken into account, these additional 
measurements are consistent with the rise time presented 
here. 
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5.2. Yet Another Second Parameter 

When spectroscopic features are taken into considera- 
tion (Benetti et al. 2005), normal Type la SNe are no 
longer believed to be characterized broadly by a sin- 
gle "first parameter" such as Amis, A, or stretch (see 
§3.2). More narrowly, a single parameter does not suf- 
fice even if attention is confined to the B and V light 
curves, which have relatively simple shapes. For exam- 
ple, Jha et al. (2007) obtain a true intrinsic color scatter 
in (B — V)35, the B — V color 35 days past tBmax, equal to 
(T(B_y)int = n !049. Careful inspection of their Fig. 8(b) 
reveals that less than 10% of o-J B _ V y mt can be ascribed 

to systematic dependence of (B — 1^)35 on the MLCS2k2 
parameter A. Therefore (B — VQ35 is a second SN pa- 
rameter. 

Broadening the discussion to include the R and I 
bands, consider as a second example the ratio TZ^ = 
isec/^s, where t sec is the interval between tg max and the 
second /-band peak, and sb is the B-band stretch. Tak- 
ing at face value the uncertainties in t scc shown in Fig. 4 
of Nobili et al. (2005), TZ^ exhibits substantial intrin- 
sic scatter while remaining largely uncorrelated with sb- 
Therefore TZ^ is another second parameter. Many more 
second parameters may await identification. 

Clearly the rise-time — fall-time difference t T { shown 
in Fig. 3, exhibiting a wl-day intrinsic scatter and no 
significant correlation with Amis, by the same criteria is 
yet another second parameter. 

5.3. Correlations with Other Variables 

To understand better the factors controlling SN rise 
and fall time, it is interesting to measure the correla- 
tions of t T , tf, and t T { with other SN variables. For 
each parameter pair, this is done here by evaluating the 
Pearson correlation coefficient R, whose Fisher function 
z' = (ln(l+i?) - ln(l-i?))/2 is distributed normally 
with variance 1/5 for this eight-SN sample. 

Mannucci et al. (2007) describe two classes of Type 
la SN bimodality that are manifested in the properties 
of their host galaxies. Used here is a simple measure 
of host-galaxy morphology: the integer T, ranging from 
—3 (E) to 7 (Sd), as shown in Fig. 6. With respect 
to T, the rise time has a positive Pearson coefficient 
i?(i r , T) = 0.671q32, while the fall-time correlation is not 

as significant, R(t{, T) = 0.46^ ; 41 . As expected, the cor- 
relation of stretch with T is intermediate between that 
of t T and tf. Sullivan et al. (2006) already have estab- 
lished at a higher statistical level the positive correlation 
of stretch with the rate of star formation. 

Also examined here are the color indices B — V, V—R, 
R—I, and linear combinations thereof. (Colors are mea- 
sured at ismax after correcting for Galactic dust using 
the map of Schlegel et al. (1998).) Studied in addition is 
vy, the initial rate of change of B — V color (§2.3). For 
this small sample, no statistically significant correlations 
of these variables with t r , tf, or t T f are evident. 

Not measured are correlations with absolute SN mag- 
nitude M: for only three of these SNe are distances well 
known (SN 2004eo is barely in the Hubble flow, and 
SN 1990N and SN 1998aq are HST Cepheid calibrated). 
Despite their 2-day difference in t rf , SN 1990N and SN 
1998aq differ in M v only by m 09 ± m 17 (Riess et al. 
2005). 
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Fig. 6. — Fitted rise time (filled circles) or fall time (open 
squares), as in Fig. 2, vs. host galaxy morphology. Both statis- 
tical and systematic uncertainties are included in the error bars. 
The positive Pearson correlation coefficient R is more significant 
for rise time than for fall time. 

Naively, a shorter rise time would be associated with a 
larger photospheric expansion velocity u CX p. Near is ma x, 
values of w cxp estimated from the Si II 6355-A minimum 
are provided for six of these SNe by Hachinger et al. 
(2006). Of these six SNe, all but SN 2002bo belong 
to the low-expansion-velocity-gradient (lvg) group de- 
scribed by Benetti et al. (2005); SN 2002bo, with a higher 
velocity gradient, belongs to the hvg group. As for the 
remaining two SNe, v cxp for SN 1998aq was observed by 
Vinko et al. (1999). Its relative stability near ts ma x sug- 
gests that SN 1998aq should also be classified as lvg. SN 
2005cf was assigned to the lvg group and its v cxp was 
measured by Garavini et al. (2007). These eight values 
of Vcxp exhibit a Pearson coefficient of correlation with 

'0.12t°ii, but the cor- 



fall time that is negligible, R 
relation with rise time is more evident, R = — 0.39^'^. 
The latter value differs from zero with the expected sign, 
but only at the one-standard-deviation level. 

6. IMPLICATIONS FOR PRECISION COSMOLOGY 

6.1. Rise Time Evolution 

For use in measuring the expansion history of the Uni- 
verse, the key aspect of any SN property is its possible 
evolution between high and low z. As pointed out by 
Nugent (1998) and Riess et al. (1999a), evidence for sub- 
stantial rise-time evolution could challenge the assump- 
tion that, after appropriate correction for light-curve 
widths, SN la absolute luminosities do not significantly 
evolve (for a recent discussion, see (Howell et al. 2007)). 
The measurement of rise-time evolution was tackled by 
Con06, who used recent SNLS data to achieve a large 
reduction in high-z rise-time uncertainty. Equally im- 
portantly for the high-z vs. low-z comparison, Con06 
applied similar analysis tools, including SN light-curve 
templates, to both samples. 
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Because the quality of high-z SN photometry, despite 
recent advances, does not yet allow accurate definition of 
multicolor light curves for individual SNe without impos- 
ing external templates, the analysis reported here does 
not compete with Con06 in comparing high-z to low- 2 
rise times. The AQUAA contribution is only to suggest 
that intrinsic variance, especially at low z, should play a 
larger role in their error analysis. 

6.2. Correcting for Light-Curve Width 

Apart from possible evolutionary effects, rise-time val- 
ues influence measurements of distance moduli through 
their contribution to the light-curve widths for which SN 
absolute luminosities M are corrected. To explore this 
point, it is convenient to divide the B light curve into 
pre-tsmax (pre-max) and post-max segments, using the 
rise time t T and fall time tf, respectively, as measures of 
the pre-max and post-max segment widths. The com- 
bination of t T and tf that effectively is used to correct 
Mb depends on the choice of first parameter (§3.2). For 
example, Amis depends only on tf, while, to a first ap- 
proximation, stretch sb is proportional to a linear com- 
bination t r + -ytf (§2.2). The value of 7 depends on the 
distribution of photometric errors over the time interval 
where the stretch fit is made; for the following discus- 
sion a representative value 7 = 2 is adopted. In §4.2 a 
« 1-day rms variation a™* in the rise-time — fall-time 
difference t T f was reported, but the dependence on t T f of 
Mb is not yet known. 

Choosing tf and t T f as the independent variables (§4.1), 
one may rewrite 



sb oc tf + 



1 



-t r 



7+1 



(11) 



For correcting the absolute SN magnitude M B , two sim- 
ple possibilities are illustrative: (i) t r and t r f are equally 
correlated with Mb] or (ii) t r f is uncorrelated with Mb- 
In case (i), correcting Mb with any linear combination of 
rise time and fall time will yield essentially the same re- 
sult. In case (ii), results from correcting Mb with stretch 
vs. correcting with fall time (which is tightly coupled to 
Am.15) need not agree. After proper calibration, these 
two choices will yield an rms intrinsic mutual difference 
in M: 



int 



int 
<7 rf 



m 033 . 



(12) 

7+1 15 days ' V ' 

where a « 1.5 is a typical stretch-correction coefficient. 
Without more information, one cannot determine which 
of these two correction methods is better. At present 
this choice is not of great practical consequence since, in 
quadrature with the « n ?18 Hubble-line dispersion, aff 
is negligible. 

6.3. Using a Template with a Biased Rise Time 

In analogy with equation (12), the systematic error 
made by correcting Mb with a stretch that is measured 
using a template having a rise time that is too long by 
+2 days is 



cf ys 



sr 



-0 m 067, 



(13) 



7+1 15 days 

which, representing w 1/3 of the magnitude of the dark 
energy signal at z ~ 0.5, is not negligible. Immunity to 



this source of systematic error often is sought by relying 
on the same template to fit high-z and low-z light curves, 
so that fitted values of high-z and low-z stretch would be 
biased in the same way. 

Unfortunately, if the high-z and low-z light curves are 
sampled differently, the cancellation in stretch bias is in- 
complete. At high z, SNe often are discovered using a 
"rolling" strategy, in which the same sky area is imaged 
repeatedly with a rapid (few-day) cadence. After a SN is 
discovered, its fluxes measured from stored pre-discovery 
images are added to the dataset, providing a uniformly 
sampled light curve. At low z, a rolling strategy is obser- 
vationally expensive due to the much larger sky area that 
must be imaged to monitor a comparable search volume; 
many nearby literature SNe were discovered and mea- 
sured in other ways. For example, of the 61 nearby SNe 
fitted by Wang, Strovink et al. (2006), 22 (9) included 
no B measurements before t s max (t s max + 5 days) . 

Consider a model experiment in which 100% of high- 
z and 60% of low-z light curves are sampled uniformly, 
while the remaining 40% of low-z light curves are sam- 
pled only after tB max - From equation (13), using a tem- 
plate with a 2-day bias in rise time would introduce a 
systematic bias in the difference of high-z and low-z ab- 
solute SN magnitudes of w n ?027. Its order of magnitude 
is that of the largest single systematic error source con- 
sidered in the SNLS first-year paper (Astier et al. 2006). 

7. PROSPECTS 

In the short term, high priority is given to build- 
ing a second, statistically independent version of the 
(presently unusual) distribution in Figs. 3 and 4. Ide- 
ally, the SNe to be fitted for this purpose would share 
features with the best in the presently available sample: 
well-observed, with BV photometry available «2 m be- 
fore -B max is reached, and sampled with < 4-day cadence 
around fg max ; not too distant (z < 0.02), to minimize K- 
corrections; not too dusty (Ay < 1), to minimize dered- 
dening corrections; and not peculiar. Full results of fits 
to shared photometry are shared with the provider. 

It would be advantageous to have available a SN light- 
curve-fitting tool that (like AQUAA) is smooth and flexi- 
ble enough to provide a statistically acceptable descrip- 
tion of unusually well-observed photometric data, but 
that (unlike the present version of AQUAA) is also able 
gracefully and objectively to represent much coarser, 
sparser data. A single tool could address these goals 
if its many fitted parameters were orthogonalized and 
mostly constrained by gentle priors derived from fits to 
the best-observed SNe. 
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APPENDIX 

A. SUPERNOVA FIT DETAILS 

This Appendix collects miscellaneous facts and observations about photometric data for the SNe analyzed here, and 
about the AQUAA fits to those data. This information is of possible interest to experts on these particular SNe. 

A.l. SN 1994-D 

The photometry used here is from Richmond et al. (1995) and Patat et al. (1996). In both papers, guidance on 
photometric uncertainties is provided, but point-by-point errors are not tabulated. Therefore the iterated fitting 
procedure (§2.5) used elsewhere for perturbing the errors was used here to set the error scale for each telescope and 
band. Within a given telescope's dataset, errors also were varied in proportion to flux -1 / 2 and in concert with the 
recorded seeing conditions. It was checked that the fitted light-curve residuals, after normalization to these errors, 
exhibit no strong flux dependence. 

A.2. SN 1998aq 

The earliest B datum is only « 1™2 below B max . Fortunately, Rie99b transformed three unfiltered additional points 
to the V band, including one that is w 2™8 below V mlix . (Not independently, Rie99b transformed these same points 
also to the B band; they are not used here, as the V filter more closely matches the unfiltered passband.) Because two 
of the transformed V points overlap well with two filtered V measurements, AQUAA was able to allow the zero point 
of the unfiltered set to float freely; the fit is sensitive only to their relative magnitudes. Nevertheless, the fitted zero 
point of the unfiltered data is in excellent agreement with that assigned by Rie99b. 

A.3. SN 1990N 

As they did for SN 1998aq, Rie99b transformed eight unfiltered additional SN 1990N points to the V band. Of the 
0™13 error they assigned to those points, 0™05 was treated in the AQUAA fit as a point-to-point error. Again, as for 
SN 1998aq, the transformed unfiltered points overlap well with two filtered V measurements, so their zero point was 
left as a free parameter in the AQUAA fit. Again it is in excellent agreement with the Ric99b value. 

Gaps in observation of ~8, ~7, and «6 days occur between iemax — 7 and tsmax + 22 days. The AQUAA fit became 
more robust when the stretch-corrected B and V spline knot phases were fixed to the values fitted for SN 2001el. The 
quoted uncertainties include an allowance for the effect of applying this extra constraint. 

A.4. SN 2001el 

The host-galaxy-absorption part of the i?-correction to SN 2001el used the values Ay = 0™54 and R v = 2.15 
measured by Krisciunas et al. (2007), who compared SN 2001cl to its less extinguished "clone" SN 2004S. 

A.5. SN 2002bo 

The host-galaxy-absorption part of the i?-correction to SN 2002bo used the value Rv = 2.6 (intermediate between 
Galactic and SN 2001el values) and a color excess E(B — V) = 0"!44 derived from the AQUAA fit. This yields A v = 1.14, 
essentially the same value obtained by the optical photometry of Krisciunas et al. (2004). 

A.6. SN 2003du 

The earliest point from HCT (Anupama et al. 2005), w 1™75 below £? ma x, lies 7.6 quoted standard deviations from 
a smooth curve drawn through nearby points from HCT and from the NOT and Asiago 1.8m telescopes (Stanishev 
et al. 2007); the other 42 points lie within 2.3cr of the fitted curve. Therefore its uncertainty was increased to ±0™15 
(see Fig. 1). In this isolated case, the assumption that the fitted B light curve is smooth led to the use of information 
from later points to reduce the relative weight of an earlier point. 

A.7. SN 2004eo 

The photometry is from Pastorello et al. (2007a) and from the preliminary B and V points shown in Fig. 5 of 
Hamuy et al. (2006). Independent fits to each dataset were performed; all quoted parameters are unweighted averages 
of the two fit values. The quoted uncertainties include allowances for possible systematic differences between the two 
components. 

The K- and .^-corrections to the data of Pastorello et al. (2007a) are taken from their paper, while those to the data 
of Hamuy et al. (2006) are discussed in §3.1. To the latter points, which are those shown in Fig. 1, uncertainties of 
0™015 were assigned at peak; elsewhere they grew in proportion to flux -1 / 2 . Again it was checked that the normalized 
light-curve residuals exhibit no strong flux dependence. 

A.8. SN 2005c] 

The B magnitude at JD 2,453,542.53 from Pastorello et al. (2007b) lies 7.1 quoted standard deviations from the 
fitted curve; the other 38 points lie within 1.8<r of the curve. Therefore its uncertainty was increased to ±0^15 (see 
Fig. 1). 
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TABLE 2 

Comparison of published B-band times of maximum flux and rise times to aquaa values 

for individual sne. 



CISJ 

DIN 


r^UDllSIlCa 

l Bmax 

(JD-244000) 


Published a 
— AQUAA 
^Brnax 

(days) 


Published 
B rise time 
(rest-frame 
days) 


Published 

A ATT A A 
— ACjUAA 

(rest-frame 
days) 


Rise time 
reference 


SN 1990N 


8083.0 ± 1.0 


0.5 


18.5 ± 1.5 b 


-1.5 


(1) 


SN 1990N 


8083.4 C 


0.9 


21.4 ±0.3 


1.4 


(2) 


SN 1994D 


9432.9 ± 0.2 


0.6 


17.6±0.5 d 


2.2 


(3) 


SN 1994D 


9432.7 C 


0.4 


15.5 ±0.5 


0.1 


(2) 


SN 2002bo 


12356.5 ± 0.5 


0.1 


17.9 ±0.5 


1.8 


(4) 


SN 2004co 


13279.2 ± 0.5 


0.5 


17.7 ±0.6 


1.1 


(5) 


SN 2005cf 


13534.0 ± 0.3 


0.2 


18.6 ±0.4 


2.0 


(6) 


Average 1 




0.44 ±0.15 




1.28 ± 0.33 




RMS f 




0.33 




0.75 





References. — (1) Leibundgut et al. (1991); (2) Riess et al. (1999b); (3) Vacca & Leibundgut 
(1996); (4) Benetti et al. (2004); (5) Pastorello et al. (2007a); (6) Pastorello et al. (2007b). 
a Value quoted by or deduced from information in rise time reference. 
b A range of 17-20 days is quoted in reference (1). 

c Estimate deduced by identifying the last two entries for SN 1990N in Table 4 of reference (2) 

with the entries for JD - 2440000 = 8071.5 and 8072.5 in Table 5 of Lira et al. (1998). 

d In a more recent analysis using methods similar to those of reference (3), Contardo et al. (2000) 

obtained a best fit bolometric rise time for SN 1994D that was 1.4 days shorter than that of 

reference (3). If the B rise time were similarly reduced, it would fall within 1 day of the AQUAA 

value. 

e Estimate deduced by identifying the first and third entry for SN 1994D in Table 4 of reference 
(2) with the first two entries in Table 5 of Richmond et al. (1995). 
Excludes references (1) and (3), in which a unique rise time function was used. 

B. COMPARISON OF PUBLISHED RISE TIMES OF INDIVIDUAL SNe TO AQUAA VALUES 

Table 2 collects published rise times for five SNe studied in this paper. Also tabulated are the times of maximum B 
flux from the published rise time analyses. Excluding the measurements of SN 1990N by Leibundgut et al. (1991) and 
of SN 1994D by Vacca & Leibundgut (1996), in which a unique rise-time function was used, on average the published 
rise times exceed those of AQUAA by 1.28 ± 0.33 days. Of this excess, an average of 0.44 ± 0.15 days may be attributed 
to the later times of maximum B flux used in the published analyses; this point was discussed in §3.3. There the 
remaining w 0.8-day excess was ascribed to two properties of the piecewise analysis method of Rie99b used in those 
publications: at ij ; n w 10 days before tsmax, the two fitted pieces of the light curve are not required to join; and, for 
t < ijoin , B — V color is assumed not to vary. Here the implications of these two properties are illustrated. 

Consider the plot of £>-band root flux vs. time shown for SN 2004eo in Fig. 1; suppose that a linear fit is made to 
early points. If, as for AQUAA, fluxes only within the first 40% of the rise time are assumed to rise parabolically, only 
the earliest three points belong in this fit. On the other hand, if all five points up to i_B max — 8.9 days are included, \ 2 
is still acceptable (4.7/3 dof), but the fitted time of explosion is earlier by 0.4 days. If they are fitted in isolation, all 
five points are not badly represented by a straight line; however, if (as by aquaa) they are fitted as part of a full light 
curve with continuous th through 3 rd derivatives, the portion of the best-fitted curve that approximates a straight 
line is shorter, as is the fitted rise time itself. 

Secondly, consider the V-band plot for SN 1990N in Fig. 1, and suppose that a linear fit is made to the points with 
r < —10 days. However, following equation (4) and the discussion in §2.1, if the B-band rise is parabolic and if the 
best-fit initial rate of change oi B — V color is negative, as it is for SN 1990N, the curve that instead should be fitted 
to these points is slightly concave downward, as shown in Fig. 1. Both fitted values of % 2 are equally acceptable, but 
the straight line yields a fitted time of explosion that is 0.3 days earlier than that of the curve. 
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